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A theory of the frequency spectrum of radio noise at extremely low frequency (ELF) is 
presented and the results compared with recent measurements of the first five "Schumann" 
resonant modes (between 8 and 34 c/s) made by Balser and Wagner [I960]. The source of 
this noise is assumed to be return strokes in vertical cloud-ground lightning flashes distributed 
randomly in time, uniformly in angular displacement along the earth relative to the observer 
and with statistics of stroke duration, interstroke intervals and strokes per flash taken from 
studies of thunderstorms reported by J. C. Williams. Thus, the mathematical model for the 
noise sources is an extremely simple one, being analogous to the shot effect in electron devices. 
The electromagnetic model employs the familiar waveguide mode theory, assumes a sharply 
bounded homogeneous ionosphere, and neglects the earth's magnetic field. 

Agreement between the shape of the theoretical and observed spectrum is good for the 
first three modes and rather poor for the higher modes. 

It is found by matching the theoretical resonant frequencies to the observed resonances 
that the product of effective ionosphere height h and the square root of effective conductivity 
V<r t - is a decreasing function of frequency. The functional dependence of this quantity on 
frequency is determined and used in the calculation of the mode spectrum. 

Discrepancies between the theory and experimental results are believed to be partially 
due to the artificiality of the sharply bounded homogeneous ionosphere model and to failure 
to give sufficient probability weighting to equatorial regions of abnormally high thunder- 
storm activity. 

These last items are the subjects of continuing work on the extension of the theory. 

1. Introduction 

The theory of extremely low frequency (ELF) propagation, as originally developed by 
Schumann [1957] and extended by Wait [1960a] and others, predicts the existence of a series of 
modes whose resonant frequencies, relative amplitudes, and q's can be found by considering the 
earth-air-ionosphere cavity as a linear filter and calculating its frequency response function 
directly from electromagnetic theory. Recently, Balser and Wagner [1960] have observed these 
modes and experimentally determined the detailed shape of the mode spectrum. 

The basis of the cavity analog is the classical Sommerfeld theory of propagation around a 
spherical earth, which leads to the representation of the field between earth and ionosphere as 
a series of "waveguide modes. " If the ionosphere height is extremely small compared 
to wavelength, as is true at ELF, the electric field is essentially radial and a single waveguide 
mode is sufficient to describe the fields at a distance from a localized source. This waveguide 
mode when considered in the time domain can, in turn, be shown (by an appropriate expansion 
of the hypergeometric series into a series of Legendre polynomials) to be a superposition of 

waveforms each closely approximating a damped sinusoid of frequency - — sJn(n-\-l), where c 

and a are velocity of light and earth radius, respectively and n is a mode index number running 
from zero to infinity. This phenomenon is entirely analogous to the ringing of a resonant 
cavity in response to impulse excitation. 

Experimentally observed VLF noise is thought to consist of a superposition of responses 
from lightning flashes all over the world [Williams, 1959]. Therefore, assuming that this 
hypothesis is correct down to ELF, a theoretical model based on it should yi eld accurate 



i Work supported by the Office of Naval Research. 

2 A preliminary note on this work appeared in the J. Geophys. Research, May 1961. 
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calculations of the spectrum of observed ELF noise. It is the purpose of the study reported 
here to construct such a model and use it to calculate the expected ELF noise spectrum. The 
results are compared with the measurements of Balser and Wagner [I960]. 

2. ELF Cavity Mode Theory 

Following Wait [1960a] one can express the significant (radial) component of the 
electric field at an observation point due to the harmonic excitation (at frequency o>) of a 
vertical electric dipole of moment Ids on the surface of the earth in the form 

E r (io>)=Idsf:C n (d)F n (ia>) (1) 

n=0 

where the ionosphere is assumed to be electrically uniform and sharply bounded at height k, 
and the following definitions apply: 

€ = Dielectric constant of free space 
P„(cos 0)=Legendre polynomial of order n 
0=Polar angle relative to the source 



and 



with 



F n (iw)z- 



and 

where 
and 



"ico[— o) 2 v (iw)+o) 2 n ] 



^(iu)^ v(v+l) = -(ia>) 2 -a(ia>) 3/2 
a 



o>l=-^n(n+l) 



Mn<Ti \ V <r g J 
Vg=(Tg-\riegUg~<rg=comi)\ex conductivity of the earth 



e 2 f N *\ 
<Ti= -^ (n _z \ X=<Ti+i«€i« (r 4 = complex conductivity 



"m\(flf— iw) 



of the ionosphere, where N, e, m, g, are the electron density, charge, mass, and collision frequency 
[Stratton, 1941], and it is assumed that co << g. In the calculation of the appropriate iono- 
spheric electrical parameters the effects of the earth's magnetic field and the motions of the 
heavy ions are neglected; in addition, the electron collision frequency is taken to be independent 
of the nonrandom electron energy. In general l^l^^l^l, so in what follows a, which will be 
called the "ionospheric loss parameter," is approximately a real quantity, equal to lj{h^ r yiff i ). 

The function E r (iu>) is actually the Fourier transform of the field in response to a (vertical) 
current impulse of magnitude / at 6=0 and at £=0, i.e., it can be regarded as the characteristic 
frequency response function of the earth-air-ionosphere cavity. F n (iu) is the Fourier transform 
of the nth mode excited by an impulse of dipole current of magnitude 7. From knowledge 
of all of the F n (io)) y s one can predict the frequency response to any time varying source current 
using standard linear system theory. 

582 



3. Superposition of Fields From ELF Sources 

The radial field component is expressed here in a spherical coordinate system wherein the 
source is at the pole. Because of azimuthal symmetry about the source, we can without 
changing the form of (1 ) consider E T (a)) as expressed in a system with pole at the observer and 
source at polar angle 6. Suppose that lightning flashes, assumed equivalent to vertical dipole 
sources, are distributed around the earth's surface. We define E r (a);d,<t>) as the Fourier trans- 
form of the radial field component at the pole due to a vertical dipole source at the point #,</>. 
The source current moment waveform at 0,<£ and its Fourier transform are denoted by i(t;d,<t>) 
and /"(a?;0,(/>), respectively. 

The voltage induced in a receiver at the observation point at time t due to the presence 
of lightning sources all over the world is proportional to 

»(*)=£ P dt'f(t-t'; e k ,4> k )j(t'; 6 k ,4> k ) (2) 

/r=lJ-co 

where the index k is used to number the source positions, from 1 through N, and 

/(«; *,*)=«S (2n+l)P n (cos e)j n (t) 

71 = 

where f n (t) is the inverse Fourier transform of F n (i cj), i.e., the impulse response of the nth 
earth-air-ionosphere cavity node. 

4. Calculation of the Power Spectrum 

The evaluation of the power spectrum of v(t) falls within a familiar class of problems in 
random noise theory, specifically, the spectral analysis of the voltage arising from the "shot 
effect" in electron devices. This problem has been treated by Rice [Wax, 1954] and others 
[Davenport and Root, 1958]. The analysis to follow is mathematically equivalent to that of 
the shot effect. 

Consider a voltage, v(t), appearing at the output of a linear filter with impulse response 
j(t) due to a train of input pulses, each of which is denoted by %(t—t k ;P). The shape of the 
kth pulse is determined by the components of a vector 0, its magnitude (time integral of ampli- 
tude over duration of pulse) is a k , and its time of occurence is t k . The output voltage is thus 

t>(*)=-S «* P Kt-t')Z(t'-t k ) fo)dt' (3) 

k U — CO 

We assume a time interval 3 — T<t< T over which all the action takes place and that there 
are exactly N pulses within that interval. Because of finite filter "memory" and the filter 
realizability requirement, 4 we can regard as zero the response to pulses occurring at values t k 
outside the interval (t— T ) to t. Here T is the time of occurence of the earliest imput influ- 
encing the filter output at time t. The interval (t—T ) is less than T but somewhat greater 
than the filter time constant. Numbering from 1 to N the pulses within the interval — T to T, 
and using (3) we can write 

*>(*)=2>*ff (*-**; ft) (4) 

fc=i 

where 

9(t-t k ) &)= P f(t-t')!;(t'-t k , $ k )dt'] -T<t k <T) fc=l, ...,N, 

Jt-T Q 

the response of the filter at time Uo a pulse occurring at time t k . 



» The time T washes out in the final results; the finite time assumption is strictly for mathematical convenience in the derivation. 
4 i.e., the requirement that the filter can respond only to past (and not to future) inputs. 
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We now assume that: (1) The a k s are independent and identically distributed random 
variables with common probability density function p(a); (2) for a fixed number of pulses 
within the interval t—T to t, the t k s are statistically independent random variables, each 
uniformly distributed throughout the interval — T to T; (3) the probability of finding exactly 
N pulses within the interval t— T to t, denoted by p(N), is independent of the statistics of the 
a k s and fi k s; (4) the fttjs are independent and identically distributed random variables with 
common probability density function p(P); and (5) the random variables a k , t k , and fi k are 
all mutually statistically independent. 

By calculating the autocorrelation function of v(t) under the assumptions (1) through (5) 
and taking its Fourier transform, we obtain for the spectrum 5 

1 



S(«)=^{W^-l)Xa>^(0>< ff (*+T)>5(«)}+iS(co) (5) 

where 

&± jAll/8-space — " - 

and 

roo 
dte-^ l g(t;0) 

\) — oo 

By the well-known convolution theorem for Fourier transforms 

fl r («;P)=F(*«)f(«;@ (6) 

where F(io>) is the frequency response function of the filter and |(co;j8) is the Fourier transform 
of the pulse. 

The first term of (5) is a "DC spike," and is not required to determine the shape of the 
mode spectrum in the frequency region of interest. We will, therefore, concern ourselves only 
with S(co). 

5. Statistics of Lightning Source Currents 

The result (5) can be adapted to the calculation of the ELF noise spectrum by allowing 
the pulses %(t—t k ;@ k ) to represent lightning flashes whose shape parameters, time of occurrence t k , 
and polar angles 6 k are assumed to be random variables, statistically independent when j^k. 
The angles B k are included among the components of the vector & k , as are the number of strokes 
in a flash [Williams, chapter VI], and the shaping parameters associated with individual strokes. 

The statistics of the lightning flashes are assumed to obey all of the assumptions (1) 
through (5). It does not necessarily follow that we can with physical justification define a 
spatial probability density function p{6) apart from consideration of other random variables 
associated with flashes. However, this assumption is basic to what follows and we note that 
it is tantamount to the assumption that statistics of lightning flash shape (i.e., those types of 
flash important as sources of VLF noise) are independent of those of terrestrial location of the 
flashes. The essential correctness of the latter assertion is supported by experimental evidence 
[Williams, tables 6-1, 6-2, 6-4, figures 6-1, 6-2, 6-7, 6-11, 6-14]. 

Using (5) and (6), the definition of f(t) given in (2), and the assumptions about lightning 
flashes discussed above, we obtain for the ELF noise spectrum (exclusive of the DO component) 

&(«)=x<iS(«,/s')i> i: flu») (7) 

— m,n=0 

where K is a constant, B' represents all lightning flash parameters exclusive of 0, and 

SU^) = (2m+l)(2n+l)F m (io:)Ft(i^L mnj 

s Ensemble averaging over /S* is denoted by <_.>, averaging over the time variables tk by a line over the symbol. 
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where 

L mn =27r I dd sin $P m (cos 0)P w (eos 6)p(6), 

(lf(co;^0t 2 )=J^(^)li(o>;^)l 2 

All p' space 

p(J}') is the probability density of the set of parameters #', and p(6) sin Odd is the probability 
that the lightning; flash is located between 6 and 6+dd. 

We are interested only in relative and not absolute amplitudes; therefore, in numerical 
computations we will fix K by assuming a normalization that equates the average power levels 
in the theoretical and experimental results. 

6. Flash Shape Parameters 

J. C. Williams [ch. V, sec. 5f] has concluded from his studies of thunderstorms that the 
overwhelming cause of VLF noise is the presence of vertical currents in the return strokes 
associated with cloud-to-ground lightning flashes. These flashes consist of from 2 to 12 strokes 
of similar shape. He has compiled statistics on flash parameters, including the shape of the 
median return stroke [ch. VI, sec. 1, 4, 5, and ch. VII, sec. 1], interstroke interval, and number 
of strokes per flash. He has also made spectral calculations for the VLF noise resulting from 
the flashes, but these are in a higher frequency range and inapplicable to our problem (ch. VII, 
sec. 3). 

The current moment of Williams' median return stroke is proportional to 

+io 3 6- ( '-w (5x ^ (8) 

where v(t) is the velocity of the lightning charge along the cloud-ground channel, the unit step 
is inserted to assure that the stroke has zero value before time t k and the lightning current 
i(t—t Jc ), v(t) and t k are measured in amperes, kilometers/second and seconds, respectively. 
The normalized spectrum of (8) is proportional to 



|/(-)! 2 = 



Mt V^ T, r A ^ \ I 



(9) 



where I t and r t are respectively the currents and time constants given in (8) and r ^=1.80(10~ 4 ) 
sec. 

It would be desirable in the interest of accuracy to consider the shaping parameters of the 
stroke model as random variables and average the spectrum over these variables. This is a 
formidable mathematical problem in general; and moreover, there is insufficient experimental 
data available to assign probability density functions to these parameters, a measure which 
would be required to evaluate the averages in the general case. The assumption that all shaping 
parameters have small standard deviations relative to their mean values would reduce the 
problem to one requiring only a knowledge of first and second moments of the distribution of 
each parameter. However, experimental knowledge of even these quantities is not available 
in enough detail to allow any meaningful assignment of numbers to the calculation of the 
spectrum of the median return stroke itself. This is tantamount to assuming that all return 
strokes are nearly identical in shape. This is not quite accurate, but the analysis will still 
provide useful insights into the effects of stroke shape on the ELF spectrum. 
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The waveform of a flash consisting of M identical strokes 6 is 

M 

*(*)=]£ *'(*-*i), (io) 

1=1 

where i(t—ti) is given by (8). 

The spectrum of a lightning flash is roughly proportional to the averaged absolute square 
of the Fourier transform of £(£), i.e., 

<ii(«)«>=i7(«)i» s (e 1 "'"-'*') (id 

where the indicated average is over the ensemble of the interstroke time intervals t t — t k . 
An obvious classification of terms in the summation of (11) results in 

"* (a) - W =M+ S • l " , - l, '<«'-'-> (12) 

where T is the average interstroke interval and ATi k is the deviation of the interval (ti—fa) 
from (l-k)T. 

To calculate (e l0}AT u) accurately, we would require a knowledge of the probability density 
function for AT ik . An extremely rough estimate of this function is provided by Williams 
[ch. VI, figs. 6-7], but it is not specified in sufficient detail to allow a precise representation. A 
rough processing of Williams' curve shows a Gaussian distribution to be a reasonable approxi- 
mation for this variable, with standard deviation a lk , the latter being assumed somewhat smaller 
than T. The result of a simple integration is 

( e ^T lk) = e -«W lk l2 (13) 

In the absence of experimental data on second order statistics of AT Uc , it will be assumed 
that the random variable (ti—(ti)) is not strongly correlated with (t k —(t k )), where k^l. 
It follows that 

<r!*~<r? forallZ^& (14) 

Invoking (13) and (14) in (12), using the power series for 1/(1— x) to bring the summation into 
closed form and averaging over all possible values of M, we arrive at the results. 

p(u) = y: p M (o>)p(M) (15) 

M=0 

where p(M) is the probability of M strokes in a flash and 



sm 2 — 



(16) 



Averaging (16) over values of M, we obtain 

P (co)-(M)(l-6-^ 2/2 )+6-^ /2 2 A n cosncoT (17) 



where 



A k = S (M-K)p(M). 

M=K 



6 Actually the strokes are not quite identical, the first being about twice as intense as succeeding ones. However, the approximation will not 
introduce serious errors into an analysis of this type. 
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The parameters T and o t have been roughly estimated from Williams (fig. 6-7) and found 
to be about 33 and 7 msec, respectively. The precision of these estimates is questionable, 
but they are probably correct to within 20 percent. The same remarks apply to similar rough 
estimates made of p(M) from Williams' cumulative probability distribution curves for the 
number of strokes per Hash (fig. 6-1). With these estimates we conclude that the probabilities 
p{\) through p{10) can be roughly approximated by the numbers 0.15, 0.22, 0.23, 0.15, 0.09, 
0.06, 0.04, 0.03, 0.02, 0.01, and 0, that (M) is about 2.5 msec, that the numbers A through 
A 9 are roughly equal to 3.5, 2.5, 1.6, 1.0, 0.61, 0.36, 0.20, 0.10, 0.04, and 0.01, and that A n 
effectively vanishes for n^>9. 

7. Geographical Distribution of Lightning Activity 

The assumption is now introduced that the probability of occurrence of a lightning flash 
within a range of angles 0, 0+d0 during the time interval between t—T and t is 



sin 



f 2r 
$p(e)dd= d4>{p(d } <f>) sin dd6) = (p +p(0))de sin 0; O<0<tt (18) 

Jo 



where p(d,(j))dd sin 6d<j> is the flash occurrence probability for a region of solid angle dd sin Bdtfr 
at 0,0, p is the part of the distribution that is uniform in the angle 0, and p(0) accounts for 
the existence of regions of abnormally high or low average thunderstorm activity. 

Initially, we will neglect p(B) and assume a completely uniform distribution p . The 
result of applying this assumption to (7) (see appendix) with the aid of the orthonormality 
of the Legendre polynomials is 

- ,,-f^, ~S (2n+l)\F n (i*)\* (19) 

Considering now the term p(Q), it has been established that there are three equatorial 
regions of abnormally high average thunderstorm frequency [Handbook of Geophysics, I960]. 
These regions are located in (1) South America, (2) Africa, and (3) the Southwest Pacific area. 
We will denote them respectively by Ri, R 2 , and R 3 . The regions center approximately at 
(50° W, -10° S), (15° E, +10° N), and (110° E, +10° N) respectively, and during all seasons 
of the year they show an average thunderstorm frequency somewhat higher than is observed 
elsewhere in the world. Speaking more quantitatively, the spatial probability density in 
those regions is from 100 to 1,000 times higher than that throughout most of the earth's surface. 

There is another near-equatorial region in Southeast Asia (^90° E, +20° N), showing 
substantial year-round thunderstorm activit}', but not as high as in R 1? R 2 , and R 3 . This 
region has not been accounted for in our calculations. Certain nonequatorial land regions 
show intense activity during their summers, e.g., the Southeastern United States (centered 
about (75° W, +30° N)), but this has also been excluded from our analysis. It is desired 
here to concentrate on those features of the ELF noise spectrum that are not season-dependent, 
and this objective is best fulfilled by studying the contributions from the three most prominent 
equatorial regions R b R 2 , and R 3 . 

In view of the approximate nature of meteorological data, it would not be feasible to attempt 
to construct a precise thunderstorm probability function for the regions under investigation. 
We will, therefore, construct an idealized theoretical model which is convenient for computation 
and which embodies the essential features of the observed distributions. 

We will approximate the probability density functions within the regions Rj hy Gaussian 
functions 

Pjid^Kj (20) 

-\Z7raj 

where Kj is a constant, 0, is the center of the distribution in R ; and aj is the standard deviation. 
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Using (7), (18), and (20), we have 

•£wm = 2^ I i ^m+^T 2_! Kjl mnj (21) 

dd sin 6P m (cos d)P n (cos 0) 



where 



v mno 



J2wa 2 j 



The integrals Z wwi have no simple orthogonality property to aid in their evaluation. The 
method to be used here will involve the expansion of the Legendre polynomials P„(cos 0) in 
a Fourier Series of the form [Jahnke-Emde, 1945] 

n 

P„(cos 0)=S ««m cos M# (22) 

where the coefficients a m are given in reference 9. From (21) and (22) we have 

m n 

where 

(0_0.)2/2cr 2 



(J mmv = dd sin 6 — j- — - — cos fid cos vd 



, 6 

V2^f 



Using trigonometric identities, assuming (justifiably, with the parameter values of regions 
Ri, R 2 , R 3 ) that (1) o-j is small compared to ir and (2) y is not near or t (so that the limits of 
the integral are effectively — oo and +°°), and integrating by completion of the square in 
the exponential, we obtain 

-4^ S (2n+l)\F„(i<S)\*=±- £ 2tK, S sin Z^-*? /2 5,(«) (24) 

2>o(|£(",/3')l 2 > n=0 F ° J=1 ' =1 

where 

B 1 („)=i („)-^ 

+ W>(«)]i + Wo(«)]i\ 

and 

4*(V)=S {[<J»+*(w)]s[«J»(«)]b 

+[&+*(«)M<Mw)] / } 

u = 

+W*-»(«)]/[A»(«)lr} 
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for &>0 where — 1, — ^— means -— 1 if & is even, — ^— if & is odd, where the subscripts R 
and /refer to "Real part" and "Imaginary part" respectively, and where 

Wi.(«)]/=S (2^+l)[F„(^)] 7 a nM 

Normalization of p(d) can be accomplished by integrating (18) subject to the above 
assumptions (1) and (2). The result is 

2^o+S Kj sin 6,e^ ,2 =*l. (25) 

j=i 

The condition (25) will serve as a constraint on the relative values of p Q and the Kj co- 
efficients. 

8. Results and Conclusions 
Note that the theoretical resonant frequencies for a lossless ionosphere are given by the 
numbers T^-JC^l »(»+l), which for c/a=46.7 are 10.5, 18.2, 25.6, 33.1, and 40.5. These 

are sufficiently different from the experimentally observed resonances to warrant the conclu- 
sion that ionospheric losses play a major role in the mode spectrum. The real positive 
quantity a, equal to 1/(AV W*) is a measure of these losses. To study its effect, the functions 
cc n F n (ico) were calculated for various values of a. In the earliest work, the mode series was 
computed for a few values of a (assumed independent of frequency) corresponding to h~80 
km (a compromise between day and night ionospheres) and values of a t between (10~ 8 ) and 
(10~ 6 )inhos/m. It was found that deviations from observed mode resonant frequencies were 
very substantial (e.g., 20 to 30 percent) and that no constant value of a would produce near- 
agreement on all five observed modes. The next approach was to allow a to vary with fre- 
quency and use the Balser and Wagner results to infer its frequency variation. This was 
initially done by graphical methods based on plots of real and imaginary parts of cx) n F n (io)) 
[Raemer, 1961.]. The method has been revised to increase accuracy. The function \o3 n F n 

a 
(ico)\ 2 was computed in terms of normalized variables x n = cx)/co n and y n —"l^- (plots are 

shown in fig. 1) and its derivative with respect to x n equated to zero to find the maximum. 
Corresponding values of y n were then found for n=l through 5, resulting in a function a(f) 
required to match theoretical and observed resonant frequencies. This function was then 
fitted to an analytical curve, by the method of least squares and (see fig. 2) a satisfactory fit 
was found to be 

a(f) = 7. 649-0. 1863/+0.01362/ 2 -0.00006602/ 3 -0.000003388f (26) 

Note that this differs from the function originally reported by the author [1961], i.e., 

«C/) = 6.44 logi /-0.0131(/-8.0) (27) 

This is due to the difference in computational method, and the function (26) should be regarded 
as the more accurate of the two. 

The fundamental assumption behind the mode-matching technique described above is 
that the observed modes are those arising from the individual nth. order frequency response 
functions of the earth-air-ionosphere cavity and are not significantly influenced either by cross- 
terms between modes, overlap between modes, or the details of the lightning flash spectrum. 
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Figure 1. \o) n F n (io3 n x n )\ 2 versus x n with y n as parameter. 
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Figure 2. Ionospheric loss parameter versus frequency. 
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Figure 3. Idealized ELF noise spectrum. 



For the moment we confine attention to the uniform angular distribution of sources. 
Using (26) to represent a, the quantity given by (19) has been calculated and is shown plotted 
in figure 3. 

The plot of figure 3 is not appropriate for comparison with the experimental curve, since 
it assumes a white lightning flash spectrum and neglects the fact (see (1)) that the field ampli- 
tude contains a factor proportional to the reciprocal of the effective ionosphere height h. If 
it is assumed that the frequency dependence of a is due entirely to the effective ionosphere 
height and not to frequency dependence in the effective conductivity <r f , then for consistency 
the power should contain a factor [a(f)] 2 . Multiplication of the quantity plotted in figure 3 
by this factor leads to degradation rather than enhancement of agreement with experimental 
results, i.e., the envelope increases rapidly with frequency. Assuming our basic model to be 
essentially correct, one of two alternative possibilities follows: (1) The frequency dependence 
in a is entirely due to h and the spectrum of lightning flashes within the frequency range of 
interest is substantially different from a "white noise" spectrum or (2) the frequency dependence 
in a is largely due to a { , in which case the lightning spectrum still cannot be " white' ' noise 
but need not decay very rapidly with frequency within the range of interest. 

Assuming that the ELF noise is largely caused by return lightning strokes, the suspicion 
that there is a significant decay in the spectrum with increasing frequency is borne out by figure 
4. This figure, showing the estimated spectrum of a lightning flash compounded from (9) and 
(17) shows an average spectral decay of between 3 and 4 db per octave of frequency in the 
region of interest. Inserting this factor into (19) results in the theoretical curve (a) of figure 5. 
Incorporating the frequency dependence of h into the calculation results in the theoretical 
curve of figure 6. On both figures 5 and 6, the experimental curve is also shown for purposes of 
comparison. In both figures the experimental curve has been normalized in such a way that 
its average power is equal to the average power of the theoretical curve. 
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Figure 4. Estimated lightning flash spectrum. 
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Figure 5. ELF noise spectrum — effective ionospheric height 
assumed constant. 
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Figure 6. ELF noise spectrum — effective ionospheric height 
assumed frequency dependent. 

The minima in the experimental curves are somewhat deeper than those of the theo- 
retical spectra and the envelope of the minima, which is roughly constant on the experi- 
mental curves, decays substantially with frequency in the theoretical results. The rate of 
decay of the envelope of the maxima is comparable, but too rapid at the high frequency 
end of the region of interest. Widths of the modes are roughly comparable. Generally the 
agreement of spectral shape is seen to be reasonably good near the lower order modes and 
rather poor near the higher order modes. On the theoretical curves of both figures 5 and 6, 
the fourth mode becomes barely visible and the fifth is completely suppressed. Agreement is 
somewhat better in general in the case where the frequency dependence of 1 /h 2 is not accounted 
for, providing some evidence that the frequency dependence of a is predominantly due to the 
conductivity <n. 

The computations on the equatorial region analysis are currently in progress. Preliminary 
results seem to indicate that accounting for these regions has no effect on the depths of the 
modes, but will have some influence on the envelope of the spectrum. 

Other extensions which might reduce the discrepancies between the theory presented here 
and the Balser and Wagner results are as follows: 

(1) The lightning statistics used here were obtained from data taken at many times of 
year and at many locations. Long term (e.g., year-round) averages of lightning parameters 
are not necessarily the same parameters averaged over a single day. The ELF noise results 
obtained by Balser and Wagner [1960] were obtained on a single day and at a single location. 
Only observed spectra averaged over many locations and times of year would be strictly com- 
parable with the theory discussed here, unless the ELF noise spectrum is found to be essentially 
independent of time and location, a supposition that would require many more experiments 
to ascertain. 
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(2) The same remarks apply to the assumption that the lightning which gives rise to 
ELF noise consists entirely of cloud-ground return strokes and that the lightning current is 
entirely vertical. Wait [1960a] has pointed out that there may be significant contributions 
from horizontal current components at the frequencies of interest here. 

(3) Assuming the return lightning stroke moment to be a random function and averaging 
over all random variables would probably make the lightning stroke more "noisy," i.e., give it 
a flatter spectrum. This could result in a theoretical spectrum whose envelope bears a closer 
resemblance to that of the observed spectrum. 

(4) The ionosphere model could be improved by accounting for the earth's magnetic field 
and the effects of heavy ions. 

(5) The probable reason for the fact that our theoretical curves are not as oscillatory as 
the experimental results is that the assumed ionospheric loss parameter is too large, tending 
to "flatten" the modes. This implies that the homogeneous ionsophere model is not sufficiently 
realistic for precise calculation of ELF noise effects. In a model accounting for the variation 
of conductivity with height near the lower edge of the ionosphere, e.g., the exponential model 
of J. R. Wait [1960b, 1960c], the loss parameter would be a complex number. 7 This would 
result in a more oscillatory spectral curve, as indicated in recent work by J. Galejs [private 
communication], as yet unpublished, and in some preliminary calculations by the author. 
Work by the author to infer information about the ionosphere from the Balser and Wagner 
curves with the assumption of a complex loss parameter is currently in progress. 



Appreciation is expressed to Dr. R. Row and Dr. J. Galejs for helpful discussions pertaining 
to the work reported here, and to Mrs. J. Van Horn and Mr. C. Roche for computational 
assistance. 

9. Appendix. Remarks on Analytical Points 

Note that with a uniform distribution, a source could be at the observer's position. The 
discontinuity of the field at this point is accounted for in the fundamental theory by the fact 
that the hyp ergeome trie series P„(— cos 0) used in the representation of the field has a discon- 
tinuity at 0=0. The expansion of this function into a series of Legendre polynomials [Magnus 
and Oberhetinger, 1949], which is basic to our formulation, is valid everywhere but at 6=0. 
Sommerfeld [1949] has shown that P v (— cos 6) becomes proportional to log 2 as 6 approaches 
zero. In carrying out the integration indicated in (7) in the vicinity of this point, the product 
(sin 6 log 6 2 ) appearing in the integrand is an indeterminate form which can be shown to vanish 
as 0— »0. The contribution to the integral at 0=0 therefore vanishes, and no difficulties are 
encountered in extending the expansion to this point. 

Another analytical point worthy of note is that the expansion mentioned above is not 
valid at u=0, ±1, ±2, ±3, . . . etc. The point u=0 occurs at zero frequency, which we 
exclude from the analysis. Also, it can be shown that, except at co=0, i; always has an imaginary 
part and therefore cannot take on values ±1, ±2, . . . etc. 
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